Intensity of parametric fluorescence pumped by ultrashort pulses 
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We investigate intensity of the parametric down conversion of ultrashort, ultraviolet pulse both 
in a low- and high-conversion regime. In the first regime, we develop a simple analytical expressions 
for the photon flux of the fluorescence. Numerical simulations using three-dimensional stochastic 
Wigner method provide results in the latter regime. We find that the perturbative approximation for 
the photon flux as a function of the wavelength of the signal and the observation angle, quantitatively 
describes even the non-perturbative regime. For short pump pulses the intensity of the fluorescence 
is highest at frequencies where the group velocities of the pump and the down-conversion are equal. 
For longer pump pulses this requirement is gradually relaxed. Additionally, for small pump beams, 
the intensity strongly depends on the spatial divergence between interacting fields. 

PACS numbers: 42.65.Lm, 42.65.Yj, 42.65. Re 
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I. INTRODUCTION 

Quantum parametric fluorescence is one of the most 
frequently used sources of nonclassical radiation. In the 
low gain regime, the parametric process is used as a ro- 
bust source of single photons and photon pairs. The 
parametric fluorescence-based sources allowed for a num- 
ber of fundamental experiments. Examples include Bell 
inequality breaking [ij, quantum cryptography [2| and 
quantum teleportation [3j. Nowadays the effort focuses 
on engineering the correlation properties [1, H, @, 0] of 
such sources. On the other hand, in some experiments 
intense pump pulses are used, which results in high para- 
metric gain. This regime is characterized by an apprecia- 
ble amount of multiple pair creation events and most con- 
veniently described by recalling the notion of squeezing 
Q . Single or two-mode squeezing described in textbooks 
can be used to describe realistic multimode parametric 
amplifier Q and it turns out that a laboratory device 
produces a superposition of squeezed states in modes de- 
fined both by the nonlinear medium and the pump pulse. 
High-gain parametric amplifiers are nowadays also fre- 
quently used for amplification of ultrashort pulses [Iol |. 
In this application the parametric fluorescence is an im- 
portant and troublesome source of noise. 

In this paper we calculate the average angular and 
spectral intensity of a spontaneous parametric fluores- 
cence. Since parametric fluorescence has quantum origin, 
its intensity cannot be calculated using classical nonlinear 
optics. Below, we adopt both perturbative and stochas- 
tic approach for this task. With a few approximations we 
provide a simple formula for the intensity of the fluores- 
cence in single-pair generation regime. Numerical simula- 
tions show, that it qualitatively reproduces main features 
of the donwnconversion spectrum in the high gain regime 
as well. 

Our results could be used as guidelines for constructing 
brighter photon-pair sources. Other application would 



be optimization of high-gain parametric amplifiers for a 
minimal contribution of parametric fluorescence pd| . 

Let us mention that a process analogous to paramet- 
ric down-conversion is observed in physics of cold atoms. 
When a pair Bose-Einstein condensates collide, a halo 
of spontaneously scattered of atoms is observed [12, EH • 
Stochastic methods used for description of this process 
proved extremely successful [HI, EH and we adopted them 
to obtain results presented below. 

The paper is organized as follows. In Section II we 
recall some basics of the theory of the parametric down- 
conversion. In Section III, we calculate the photon flux 
in the perturbative regime. In Section IV, we compare 
these results with numerical predictions in the high-gain 
regime. We conclude in Section V. 



II. THEORY OF THE PARAMETRIC 
DOWN-CONVERSION 

A process of parametric down-conversion takes place 
in a crystal with non-zero second order nonlinear sus- 
ceptibility tensor x when it is irradiated with a pump 
beam. During pass through the nonlinear medium each 
of the pump pulses gives rise to photon pairs which form 
the parametric fluorescence. Below, we formulate a the- 
oretical model describing the nonlinear interaction which 
allows for calculation of the fluorescence intensity. 

We assume, that the medium interfaces are parallel to 
x—y plane, while fields propagate along the z-coordinate, 
as depicted in Fig. [T] We expand all the fields involved 
in the interaction using an approach from nonlinear op- 
tics. At each z an expansion in a basis of plane waves 
with frequency to and wave- vector k x , k y is given, while 
z parametrizes the development of the evolution. To 
shorten the notation, let's denote the triple (oj,k x ,k y ) 
parametrizing each plane-wave component of the fields 
by k. For concreteness, we consider Type-I interaction 
in BBO (/3-Barium Borate) in which pump propagates as 
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FIG. 1: Scheme of a measurement of the photon flux n(A, a) of 
the spontaneous parametric fluorescence. The nonlinear crys- 
tal is pumped along z-axis. Out of the emerging parametric 
fluorescence we observe a small portion at the wavelength A 
propagating at an angle a to z-axis in the x-z plane. 



a basis of plane waves parametrized by k — (uj,k x ,k y ), 
which corresponds to the following expression for the field 
operator E(x, y, z, t) in the position space: 



E(x,y,z,t) 
d 3 n [ 



huj 



(2^) 3 V 2e ^R 



Here a(«, z) form a complete set of bosonic annihilation 
operators at each z while n(uj) is the index of refraction 
of the medium. The evolution of the parametric fluores- 
cence is also largely determined by the crystal dispersion, 
which determines k z for each plane wave component with 
given k = (lo, k x , k y ) at a value k z = k z (n), 



k z (n) 



i-2 



1,2 

K y . 
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an e-ray while the fluorescence as an o-ray. Results can 
however be easily generalized for any nonlinear crystal 
by substituting its particular dispersion relation in the 
formulas below. 

Let's consider a regime in which the pump consists of 
intense pulses described by a classical envelope A p (n p , z) 
in the frequency domain, spanned by k p = (tu p , k px , k py ). 
The evolution of the pump pulse A p (n p ,z) along z is 
determined by the dispersion relation of the crystal. A 
nonlinear interaction with the fluorescence field can be 
neglected, since the intensity of the pumping field is large 
and its depletion due to parametric down conversion is 
minimal. This condition is readily satisfied in many ex- 
periments. Note that the intensity of the pump is also 
assumed large enough to neglect any quantum fluctua- 
tions and describe it by a classical field. 

For a pump propagating as an e-ray in a crystal cut 
at an angle 9 to the z-axis the dispersion relation of the 
crystal reads: 



w p [k px cos 9 — k pz sin 9)' 



nl{uj p ) 



(k pz cos 9 + k px sin 9)' 1 



k 2 

— — (1) 



where n {ui p ) and n e {u> p ) are ordinary and extraordinary 
refractive indices of the crystal. For each plane wave com- 
ponent of pump field with given k p = {uj p , k px , k py ) the 
above relation fixes the z-component of the wavevector 
k pz at a value k pz (K p ). During propagation along z each 
of the plane-wave components will acquire phase equal 
to exp[ik pz (Kp)L]. Thus, the equation of evolution of the 
pump pulse A p (k Pi z) can be written as [1 €31 ] : 

d ~ 

^— A p (np 7 z) = ik Z:P (K p )A p (Kp, z). (2) 

On the contrary, the parametric fluorescence has to be 
described by a quantum field. We expand this field in 



Additionally, the signal field interacts with the pump 
pulse via the second order nonlinearity. Classical equa- 
tions, describing propagation of the fluorescence field 
through the crystal are well known [l7l | . Quantization of 
the signal field corresponds to replacement of c-number 
amplitudes and their complex conjugates of plane-wave 
modes by annihilation and creation operators a(n, z) and 
a'(K, z) and leads to the Heisenberg equations for their 
evolution: 



—a(n, z) = ik z (n)a(K, z) 



, A p (n , z) ± , 
an a' (k —k,z). (5) 

An 



Here, A is the maximal amplitude of the pump pulse in 
time-position space Aq — J d 3 k p A p (k p ) , while Lnl is a 
nonlinear length, a constant which determines a charac- 
teristic length over which nonlinear interaction produces 
a few photon pairs, defined as: 
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(6) 



For a crystal of a length L and given pump pulses 
A p (k,z = 0) at the crystal entrance we can find Green 
functions C(k, k') and S(k 7 k') of the Eq. ([5|) and express 
the fluorescence field at the crystal exit face a ou t(/i) = 
a(n, L) in terms of that on the input face ai n (n) = a(«, 0): 



t(«) 



dn' \c(k, kOMk') + S ( K > "OoLOOl ■ ( 7 ) 



The above equation gives a complete quantum descrip- 
tion of the parametric amplifier. In our case it is highly 
redundant, since we are interested in the average inten- 
sity of the parametric fluorescence (a' i '(/t)a(K)) = (n(K)) 
as a function of the frequency and direction, both con- 
tained in k. In the following section, we use Eq. only 
as a backbone for developing more specialized results. 



III. PERTURB ATIVE APPROXIMATION 



In this section we discuss the properties of the average 
intensity of the parametric fluorescence calculated in the 
first order approximation with respect to the strength 
of the nonlinearity 1/Lnl- This approximation physi- 
cally corresponds to a situation when the fluorescence 
field consists of a single or a few pairs of photons emitted 
into distinct frequency or angular regions. It is expected 
that as soon as the pump becomes so intense that the 
photons are emitted in bunches, perturbative expansion 
fails. However, as we show, some results inferred from 
the perturbative approach, can be directly applied even 
in the high-gain regime. 

Basic expressions derived in this section are well known 
and have been introduced in works on spontaneous para- 
metric down-conversion in the context of photon sources 
[J. However, the intensity of such sources have, surpris- 
ingly, received very little theoretical attention. 

In order to calculate n(re) = (a^ (k)ci(k)) which de- 
scribes the intensity of the parametric fluorescence, we 
first obtain an approximate Green functions C(k, k') and 
S(k, k). This is done by transforming into the interaction 
picture, which is accomplished by substituting: 



z) = cxp[ik z (K,)z]di(K, z), 



(8) 



into Eq.([5|. The resulting equation of evolution of 
ai{n,z) contains only nonlinear term and can be eas- 
ily solved perturbatively with respect to 1/ Lnl [9] • The 
first order solution is next compared with a general form 
given by Eq. ([7]). The C(k,k') function in the lowest 
order describes simple dispersive evolution of the fluores- 
cence filed passing through the crystal: 



C(k, k') = 5{k - n')e lk ^ )L , 



(9) 



while the S{k, k!) describes the process of pair genera- 
tion: 



xstac (i^>), (10) 

where Ak(n, k') is a wavevector mismatch between plane- 
wave components of the fluorescence characterized by 
spatio-temporal frequencies n and n' and coupled plane 
wave component of the pump characterized by spatio- 
temporal frequency k + k! : 




FIG. 2: (color online). Observation of a plane- wave 
monochromatic component of the fluorescence field, directed 
to the upper right in the picture, leads to projection of the 
conjugate idler field into a pulsed form, directed to the lower 
right. The idler wave-packet takes up all the frequency and 
wavevector span of the pump pulse, incoming from the left in 
the picture. 



vector k: 



n{n) = (al ut (K)a ou t(n)) = / d 3 k' \S(k, k')\ 2 = 



A p (n + k') 



An 
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(12) 



The above integral can be given a meaningful physical 
interpretation. Observation of a signal photon testifies 
generation of the idler photon. In our case the first is de- 
tected using narrow-band filters and is nearly plane- wave 
and monochromatic, characterized by vector n. Thus, all 
the frequency and direction distribution present in the 
pump is transferred to the idler photon, as depicted in 
Fig. [2 Each plane-wave component of the idler charac- 
terized by k' can support the production of signal char- 
acterized by k if the pair in question is phase-matched 
with a pump component characterized by k + «/, i.e. the 
phase mismatch LAk(n, k') is small. 

The integral given in Eq. (TP2)) cannot be evaluated an- 
alytically in case of pulsed pumping, when A p (n p , z) is 
nonzero over a finite range of n. However, the integrand 
can be approximated with a Gaussian function. First, 
let us assume a Gaussian form of the pump amplitude at 
the crystal entrance face A p (k p , z = 0), corresponding to 
a pulse of duration t p propagating in a beam of width 
w p , centered around frequency 2wq: 



Ak(n, k') = k pz (n + «') - k z (n) - k z (n') (11) 

The above result can be used to calculate the intensity 
of spontaneous parametric fluorescence emitted at a cer- 
tain direction in a certain frequency, corresponding to a 



A p (k p , z = 0) = A Q 



W 2 p T p 



x exp 



(2tt) 3 / 2 




(13) 
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FIG. 3: (color online). Angle a s = &rcs'm{ck° (ui) / u>] at which 
perfectly phase-matched fluorescence at the wavelength A = 
2tic/ijj emerges outside the crystal. Plot corresponds to type-I 
phase-matching in BBO pumped at 2ttc/(2ujo) — 400 nm for 
crystals cut at angles 9 = 29° (solid line), 31.3° (dashed line), 
35° (dotted line) and 40° (dash-dot line). 



Concentration of the pump energy around given k p to- 
gether with requirement of phase matching Afe(/c, k') < 
1/L reduces the range of spatio-temporal frequencies n 
and k' over which a significant amount of parametric 
fluorescence is generated. To analyze this effect, let 
us consider perfectly phase-matched parametric down- 
conversion of the central plane-wave component of the 
pump with k p — (2wo, 0, 0). This part of the fluorescence 
consists of a plane wave component-pairs characterized 
by k and k! which satisfy the following set of equations: 



k + h' = (2cj ,0,0) 
Ak(K,K r ) = 0. 



(14) 



The first equation corresponds to a requirement that 
both components of the down-conversion couple to the 
central component of the pump, while the second equa- 
tion is the perfect phase matching condition. This is in 
fact a set of four equations for six components of k and 
k! which defines a two-dimensional surface. Assuming 
k = (u>, k x , k y ) as before, from the first equation we get 
k! = {2uj$—Lu,~k x ,—ky). Then the solution of the second 
equation can be cast into the following form: 



(15) 



This is an axially-symmetric surface in three-dimensional 
space of k = (u),k x ,k v ), describing plane- wave compo- 
nents of the parametric fluorescence which can be pro- 
duced from the central component of the pump in a per- 
fectly phase-matched process. In Fig. [3] we plot a sec- 
tion through this surface in the k x -u plane, for para- 
metric down conversion in a BBO crystal pumped at 
2ttc/ (2ojq) — 400 nm and four different cut angles 9. 

As the spatio-temporal frequency of the fluorescence k 
deviates from the perfect phase matching surface k°(u>) 
its intensity diminishes. In the region of k in which the 



fluorescence intensity is large, the sine function can be 
approximated by a Gaussian [l8j]. Such an approxima- 
tion fails only for very thick crystals, in which intricacies 
of the phase matching may play important role, but the 
fluorescence emitted in such cases contains very compli- 
cated correlations, which render it useless in most of the 
known applications. 

The next step towards Gaussian approximation of the 
integrand in Eq. (|12[) is a linearization of the phase mis- 
match Afe(/c, k') given by Eq. (jTTJ) around points on per- 
fect phase-matching surface. Since it is axially symmet- 
ric, we choose points in k x > and u) half-plane. For 
each fixed observation frequency tu h s we expand phase 
mismatch Ak(K,K') around kq — (u) o hs,k°(u)),0) and 
Kq = (2ujq — w bs, — 0) up to linear terms. This 
gives 

Ak^ n ( K , K ') = (2uj -Lj'-Lu ohs )Al3' 1 + (k x -k (Lu ohs ))Ap x 

+ kyApy + (k' x + k (iJ obs ))Ap' x + k' y A P ' y , (16) 

where the expansion coefficients equal: 
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(17) 



The above coefficients have simple physical meaning. 
A/3[ is a difference of inverse group velocities between 
the pump pulse and the fluorescence component charac- 
terized by k' , while Ap give the amount of spatial dis- 
placement per unit length between pump pulses and the 
fluorescence component characterized by k or k' in x or 
y direction. Note, that all the above coefficients depend 
on the choice of the observation frequency w b s - 

Finally, in Eq. (|12p we approximate the sine function 
by a Gaussian 



sine 



f LAk z 



exp 



~(Ak^f\ (is, 



and obtain the following expression for the fluorescence 
intensity 



h(k) 



L 

x /' dVe-^[ Afe ""( K < K ')] 2 



A p (k + k') 



A 



(19) 



This integral can be evaluated analytically. Let us sub- 
stitute k — kq — (uj, k°(uj), 0) which describes the perfect 
phase-matching line. Then we obtain an estimate of the 
intensity of the fluorescence emitted at a frequency lo at 
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FIG. 4: (color online). Square idler spatial displacement co- 
efficient Ap\, + Apy/ in BBO for four different crystal cut 
angles = 29° (solid line), 9 = 31.3° (dashed line), = 35° 
(dotted line) and = 40° (dot-dashed line). 



an angle a — arctan(cfc (u))/w) to the z axis in z 
plane: 



n(wobs,&°(wobs),0) 



(- 



Air 3 / 2 \L N l 
A P ; 2 ) / LAf3[ 



-1/2 



(20) 



Let us consider the dependence of the intensity of 
the fluorescence along the perfect phase-matching line 
^(w bs, fc (w o bs), 0) on the frequency w b s at which we 
observe it. Note again, that observation of the fluores- 
cence signal at the frequency u> ^ s testifies generation of 
the conjugate (idler) photons at the frequency 2loq — oj b s - 
It turns out, that the propagation properties of the idler 
photons influence the intensity of the signal. In partic- 
ular three parameters, which vary with w bs come into 
play. The first is the ratio of the temporal walk-off be- 
tween pump pulse and conjugate photons over the crystal 
length LA/3[ to the pump pulse duration t p . The other 
two measure the ratio of the transverse displacement be- 
tween the pump pulse and conjugate photons LAp' x and 
LAp' y to the pump pulse width w p . 

This can be understood by reconsidering picture drawn 
in Fig. [5] As we argued before, observation of a nearly 
plane-wave signal transfers the frequency and wavevec- 
tor distribution of the pump to the idler photon, which 
takes a form of a short pulse and propagates according to 
the crystal dispersion. Usually, it quickly diverges from 
the pump pulse — the expression in the square brackets 
in Eq. (|2"0"|) is large. However, for some frequencies of 
the signal w bs the idler propagates with the pump for a 
longer distance which increases the nonlinear interaction 
and the photon flux n(w bs, fc (w o bs), 0). 

We plot relevant material coefficients in BBO in case 
when it is pumped at 400nm and cut at four character- 
istic crystal cut angles 9 in Fig. 2] and [5l 
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FIG. 5: (color online). Inverse of the group velocities of pump 
(dotted line), signal (solid line) and idler (dashed line) for 
following angles of the crystal setting: (a) — 29°, (b) = 
31.3°, (c) 35° and (d) 40°. 
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We also plot the intensity of the fluorescence along the 
perfect phase-matching line n(w b s , & (w o b s ), 0) given in 
Eq. (|20|) as a function of the wavelength of observation 
27rc/<x> bs in Fig.[5]for a few typical pump pulse durations 
and widths. 

It is worthwhile to consider the features displayed by 
this plots. For short pump pulses, when temporal walk- 
off between pump and conjugate photons is more pro- 
nounced, the fluorescence intensity reaches local max- 
imum values at a wavelength where this walk-off van- 
ishes A/3i = 0. For longer pulses the temporal walk- 
off loses its importance, and the intensity of the fluores- 
cence increases with wavelength due to reduced spatial 
displacement between the pump and conjugate photons, 
i.e. Ap'i + Ap' 2 . 



IV. THREE DIMENSIONAL STOCHASTIC 
MODEL 

In the previous section we calculated the intensity of 
the parametric fluorescence in the single pair genera- 
tion limit. Now we switch to intense-pumping regime, 
when photons are generated in bunches and perturbative 
approximation breaks down. Estimation of the fluores- 
cence intensity is of a great importance in novel practi- 
cal applications of parametric amplifiers. For example, 
when using this kind of nonlinear interaction to amplify 
broadband chirped pulses, the parametric fluorescence is 
a source of large and difficult to avoid noise [ll[ . 

In this section we describe a numerical method of 
solving the dynamics of a coupled pump-signal system 
given by Eqs.Q and ([5]) and estimating the intensity of 
the parametric fluorescence. We use so-called stochas- 
tic Wigner method, which can in principle provide exact 
solutions for linear quantum systems, such as described 
by Eq. (f5|). If the Wigner function of the field at the 
entrance is Gaussian, in particular representing vacuum 
state, the Wigner function of the final state of the para- 
metric fluorescence will be Gaussian too. It is easily de- 
duced from a general form of the input-output relations 
given by Eqs.([7|). 

A crucial observation is that if the Heisenberg equa- 
tion is linear, the evolution of the Wigner function is 
the same as a classical Liouville equation for the density 
probability of finding the system in a state with partic- 
ular electric field. Thus the Wigner function VF[a(ft, z)] 
can be interpreted as a probability density of finding the 
system with a fluorescence field described by a spectral 
amplitude a(n, z). 

Therefore, to find the Wigner function of the fluores- 
cence at the crystal exit face we repeat the following 
procedure. First we draw an input spectral amplitude 
a(K, z = 0) with probability given be the vacuum Wigner 
function 



W[a(n,z = 0)] = 




n(a.u.) 




FIG. 6: (color online). Photon flux of the fluorescence 
n(oj b s , fc (a-> o bs), 0) along the perfect phase matching line plot- 
ted in Fig. [3] calculated using Gaussian approximation for 
2mm long BBO crystal. Plots are made for four different 
crystal cut angles: (a) = 29°, (b) 31.3°, (c) 35° and (d) 
40°. Various lines indicate different pump pulse parameters: 
solid lines correspond to t p =60 fs, w v =80 ^m, dashed lines 
correspond to larger pump beam, t p =60 fs, w p =160 /im and 
dotted lines correspond to longer pump pulses, t p =120 fs, 
w p =80 /im. 
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FIG. 7: Graphical illustration of the Wigner method. 



where we assumed a finite number N of discrete modes 
indexed by k. Next we solve classical equations of prop- 
agation of the signal field: 



d 

—a(K, z) = ik z (n)a(K, z) + 
oz 



(Ik 



,A p {k, z)q*(k — n' , z) 



This way we obtain a possible output spectral ampli- 
tude a(n, z = L). By repeating this procedure many 
times, we obtain a set of probable output field ampli- 
tudes {a(n,z~L)}, which approximates Wigner func- 
tion W[o.(k, z = L)]. This is schematically illustrated in 

Fig.m 

This can be used to calculate the mean photon flux at 
a given direction and frequency. It equals a mean |q;(k)| 2 
averaged with the Wigner function W[a(/e, z = L)j mi- 
nus vacuum average. Using the set of output amplitudes 
{(x(k, z = L)} this is calculated as: 



7i(k) 



(23) 



We performed the procedure outlined above for inter- 
action in 2mm long BBO crystal. The pump intensity 
was adjusted so that the total number of photons in the 
fluorescence field is near 10 8 . This corresponds to the 
typical conditions found in noncollinear-OPA (NOPA) 
used for amplification of ultrashort pulses [10( . 

For each of the various pump pulse durations, widths 
and crystal cut angles we performed ten independent sim- 
ulations of classical fluorescence evolution. The average 
photon flux was computed according to the formula (|23[) . 
Additionally, we assumed that n(n) is almost indepen- 
dent of the angle around z axis and we averaged the re- 
sults over this angle. This last step made ten independent 
simulations sufficient to calculate smooth intensities. 

Figure [5] shows the results of the simulations. The 
area of high intensity winds around the perfect phase- 
matching curve pictured in Fig[31 as we argued before. 



Moreover, qualitatively, the intensity along this curve 
changes as predicted by the perturbation theory and plot- 
ted in Fig. This can be considered an unexpected 
success of the perturbation theory and qualitatively ex- 
plained as follows. The process of parametric amplifica- 
tion with intense pumping leads to intense fluorescence 
fields, which exhibit mainly classical properties. In par- 
ticular 1/2 in Eq. (f2"3"|) can be neglected, which makes 
the process virtually identical to a classical amplification 
of a weak, white noise. In the course of interaction the 
components of the fluorescence become pairwise corre- 
lated. To some extend, the noise-like fluorescence can be 
separated into pairs of peaks which have approximately 
opposite k x and k y and frequencies summing up to the 
frequency of the pump, 2u>$. Most of such pairs quickly 
diverge away from the pump pulse, but some propagate 
with it enjoying exponential growth. This is very similar 
to the mechanism behind the intensity of the fluorescence 
in a single pair generation regime, hence the intensity 
profiles are very similar. 

It is worthwhile to compare the numerical results from 
Fig. [8] with Figures [4] and \5\ We find that for a short 
pump pulses the temporal walk-off between the pump 
and the fluorescence must be very small to achieve in- 
tense generation. For longer pulses, this requirement is 
relaxed and our calculations reveal broadband paramet- 
ric generation. 



V. CONCLUSIONS 

In this article we adopted perturbative and Wigner 
stochastic approaches to the task of calculating an inten- 
sity of the parametric fluorescence. In a single pair gen- 
eration limit, some analytical expressions for the down- 
converted field can be derived. By linearizing the phase- 
matching function we found simple formulas for the in- 
tensity of the fluorescence in this regime. Photon flux of 
the signal depends on the temporal and spatial walk-off 
between the idler and the pump. In case of type I in- 
teraction in BBO this leads to peaks in the spectrum of 
the downconversion signal at the frequencies where group 
velocity of the idler matches up with the pump. 

In an intense pumping case we performed a number 
of numerical simulations and found that the spectrum 
of the fluorescence exhibits qualitatively the same fea- 
tures as predicted by the perturbative approximation. 
The contrast between the wavelengths at which the flu- 
orescence is intense and those at which is weak is more 
apparent due to exponential character of amplification in 
non-perturbative regime. 

The above results could be used as guidelines for con- 
structing more bright photon-pair sources. They could 
be also used for optimization of parametric amplifiers, in 
which one seeks minimal contribution from the paramet- 
ric fluorescence. 
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FIG. 8: The contour plot of the averaged photon flux of the fluorescence as a function of the angle a = axcsm(c\k x \/ui) and 
the wavelength A = 2-kc/ui for a 2mm long BBO. Each row corresponds to different crystal cuts: 8 = 29° (1st row), 6 = 31.3° 
(2nd row), 8 = 35° (3rd row) and 8 = 40° (4th row). Each column corresponds to different pump parameters: column A is for 
r p =60 fs, w p —80 nm, column B is for r p —60 fs, w p =160 pm and column C is for r p =120 fs, w p =80 /nm. The contour lines 
were drawn at 0.25, 0.5 and 0.75 of the maximum of each plot. 
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